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\0 • J/ip measurements in p+p and d+Au collisions serve as crucial references to understand 

^ ■ the J/ip production in Au + Au collisions at RHIC where quark gluon plasma (QGP) is 

> . expected to be formed. They also provide important clues to study various interesting 

^ ! phenomena such as the gluon shadowing and color glass condensate. We report the latest 

Q I results from PHENIX experiment on J/ ijj production in p + p and d + Au collisions at 

O ' forward, backward and midrapidity. 

O . 1. Introduction 

X : 

<D ■ J/ip enhancement [H] or strong suppression is considered to be the signature of QGP 

! formation. To quantify these effects, one needs measurements in p+p and d+Au collisions 
^ I to understand many cold nuclear effects such as gluon shadowing, nuclear absorption, 
^ I initial state energy loss, Cronin effect and gluon saturation [ 3^ that can have interesting 

■ impact on the J/ip production. 
^ ■ PHENIX[2] measures J/ip via dielectron channel using midrapidity (|?7| < 0.35) central 
arm spectrometer and dimuon channel using forward (1.2 < rj < 2.4) and backward 
(—2.2 < 7] < —1.2) rapidity muon arm spectrometer. 

In year 2003, RHIC delieved 2.74 ub'^ d + Au and 350 nb^^ p + p collisions. PHENIX 
central (muon) arm collected about 400 (1400) J/ip in d + Au and 100 (420) J/ip in p + p 
collisions. These are the first measurements of cold nuclear effect at RHIC[^. The final 
results are reported here. 



2. Results and discussion 

Figure E^a) shows the J/ ijj rapidity distribution in p + p collisions. The dashed error 
bars represent systematic uncertainties relevant for comparing the two rapidity bins in 
each muon arm, while the solid error bars represent the overall uncertainties relevant 
for comparing points at different rapidity. A 10% of overall normalization error is not 
included. The total cross section is calculated via fitting the data with the PYTHIA 
predicted rapidity shape using GRV94HO parton distribution function (PDF) and the 
J/ip dilepton decay branch ratio of 5.9%. The uncertainty from the PDFs is found to 
be less than 3% after repeating the calculation using PDFs with very different shapes. 
The cj//'^ = 2.61 ± 0.20(fit) ± 0.26(abs) /xb, where the first error is the fitting error 

*for the full list of PHENIX authors and acknowledgements, see Appendix 'Collaboration' of this volume. 
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and the second one comes from the uncertainty of our minimum-bias trigger bias. The 
result is consistent with our measurements in year 2002 [ |n|. Figure Wljo) shows the 
minimum-bias RdAu = (^dAul^ x 197 x cjpp) versus rapidity (12% overall normalization 
error is not included). The forward (backward) rapidity corresponds to Bjorken variable 
Xftj ~ 0.003 (0.01), where one expected shadowing (antishadowing) to be effective. The 
midrapidity corresponds to ~ 0.01 where minimum nuclear effects are expected. The 
EKS98 prediction [ [TO] assuming Imb absorption cross section fits the data best. The 
models assuming large shadowing effects like FGS [02] and coherence-length model [[7] 
are disfavored. 
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Figure 1. (a): Rapidity distribution in p + p 
collisions, (b): The minimum bias i^d^i versus 
rapidity. 



Figure 2. a versus pt at different 
rapidity compared to E866 measure- 
ments. 



Figure |2l shows the nuclear modification factor a [odAu = crpp{2 x 197)°) versus px- 
An additional 0.02 overall uncertainty is not shown. The dashed lines are simple fits to 
the E866 results[IS]. One can see a similar tendency as shown in E866 results that the a 
increase towards larger px in the forward and backward region. This is consistent with the 
fact that < > in p -|- p collisions (2.51 ± 0.216*6^^/ c^) is much smaller than in. d + Au 
collisions (forward: 3.63 ± 0.25Ge^^/c^, backward:4. 28 ± 0.316*61^^/0^) as expected from 
Cronin effect. One interesting observation is that the central arm measurements show no 
broadening. The midrapidity < p]^ > in p + p collision (4.31 ± 0.85Ge^^/c^) is larger by 
almost two standard deviations than that in the forward and backward rapidity. It could 
either come from statistical fluctuation or indicate something interesting. The results 
from the high luminosity 2005 p + p run should give a clear answer on this. 

Figure Efa) plots the a against Xhj in different collision energies. If shadowing or 
antishadowing is the dominant nuclear effect, one would expect that alpha scale with Xbj 



since PDF is the function of x^ and where 



m 



for J/ip production. This 



is in clear contradiction with the data. One possible explaination is that the coherence 
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length increases with colhsion energy [ |^ . The initial state energy loss dominate the 
nuclear effect at low energy (NA3) since the coherence length is short and the parton 
has large chance to lose energy before the hard scattering happens. At RHIC energy 
shadowing is the dominant effect. We also plotted the alpha against x^? as shown in figure 
Ofb) in different collision energies. The result shows interesting scaling behavior. Some 
theoretical calculations try to explain the data with the Sudakov effect [Hj. 




Figure 3. a versus (a) xa and (b) xp with 
present 200 GeV J/ip results compared to 
lower energy results [E]- An overall uncer- 
tainty of 0.02 in our a values is not shown. 



Figure 4. RdAu versus centrality at three 
different rapidity regions, compared to 
calculations [nni including absorption and 
EKS (solid) or FGS (dashed) shadowing. 



Figure |3] shows the nuclear modification factor RdAu versus centrality represented by 
the average number of collisions in each centrality bin at different rapidity regions. The 
bars at the low end of each plot represent the scaling error in each rapidity region. An 
additional 12% global error bar is not shown. The histograms at the bottom of the lower 
panel indicate the distribution of the number of collisions for each of the four centrality 
bins. One can see that in the small Xbj region, the suppression become more significant 
towards more central collisions but the overall dependence is weak. Similar but weaker 
centrality dependence is observed at midrapidity. No clear centrality dependency is found 
in the large Xbj. Clearly the FGS overpredicts the suppression in small x^j region. 

One can deduce the baseline RauAu at forward and backward rapidity via multiplying 
RdAu results in forward and backward rapidity. Similarly RauAu at midrapidity can be 
deduced via multiplying RdAu at midrapidity by itself. The RauAu calculated this way 
as shown in figure El does not fully account for some nuclear effects like the final state 
comover suppression but it tell us the amount of shadowing and absorption expected to be 
seen in Au + Au collisions. Theoretical calculations shown in figure H^b) are also plotted 
after some rescaling to approximately describe the data. Same calculation is done to get 
the RauAu prediction from these theory curves. One can see the interesting change of 
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RauAu shapes corresponding to different amount of shadowing and antishadowing effect. 
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Figure 5. RauAu deduced from measured 
RdAu at different rapidity region. 



Figure 6. First T measurements at RHIC. 



3. Summary and outlook 

In year 2003, PHENIX made the first measurement of cold nuclear effect on J ftp pro- 
duction in ^/SNN = 200 GeV d+Au collisions at RHIC. The results show weakly increased 
J/ip suppression towards more central collision and smaller Xbj region, px broadening is 
observed although the midrapidity results show no such effect possibly due to the statsit- 
ical fluctuation. More luminosity of ci + Au collisions is needed to do detailed study on 
cold nuclear effect. In year 2005, PHENIX accumulated much higher luminosity of p + p 
collisions than that in year 2003 as being illustrated by first measurements of T at RHIC 
shown in figure IHl This will significantly reduce the uncertainty of our J/ip measurements 
in cold nuclear medium. 
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